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ABSTRACT 


We have studied the feasibility of building a biomass gasification plant with an innovative spouted bed reactor for 
distributed energy production. The process was simulated using a thermodynamic approach (concentrated param¬ 
eter model) using LIBPF, a C++ process modelling library. A nominal size of about 100 kW total thermal power was 
chosen. 

Fifteen chemical species were taken into account (dry and ash-free (daf) biomass, C, CH 4 , CO, CO 2 , H 2 , H 2 0, N 2 , 
NH 3 , 0 2 , methanol, ethylene, phenol, naphthalene and ash). 

After preliminary simulation trials and unit operation sizing, a proper plant configuration was developed. 

The air ER (equivalent ratio), syngas condensate recycling and biomass moisture fraction were studied in compli¬ 
ance with the constraint of the gas flow required for a stable spouting regime (120 m 3 /h) and optimal gasification 
temperature (930 °C). 

A working condition map was obtained in terms of air ER and syngas condensate recycling for different biomass 
moisture fractions and the main plant performance parameters were studied. 

The net electrical yield was found to vary from the ideal maximum of 25% to a minimum of about 17% and the 
heat cogeneration yield was found to vary in the range 45-48% for all the process conditions studied. 

© 2012 The Institution of Chemical Engineers. Published by Elsevier B.V. All rights reserved. 
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1. Introduction 

Today world energy consumption is about 12.5 x 10 9 toe, about 
three times that of 40 years ago (IEA, 2011), and it is still 
increasing, demonstrating that energy production is one of 
the most important demands in all industrialized countries 
(such as the those of Europe and the USA) and industrializing 
ones (such as China and India). Economic and pollution prob¬ 
lems, fossil fuel limitation and quality of life are all linked to 
the issue of energy production. A good strategy is important to 
meet the current growth in energy consumption, political sta¬ 
bility and environmental problems. In this scenario, renewable 
energy sources have a central role and the present work is ded¬ 
icated to the study of biomass as an energy source feedstock 
for an innovative thermochemical process. 

Obtaining energy from biomass in a safe way could reduce 
dependence on fossil fuels and energy importation, reduce 
damage to the environment, achieve the goal of a zero carbon 
foot print, enhance the use of energy refuse and residues, the 


return of abandoned energy sites to agricultural production 
(for energy crops), and create new business and occupational 
opportunities. 

The biomass considered in this work is wood chips, which 
are suitable for high-temperature thermochemical processes 
such as pyrolysis, gasification and combustion (Milne et al., 
1998; Higman and Burgt, 2003; Goyal et al., 2008). 

In terms of gas feed, gasification can be seen as a 
compromise between endothermic pyrolysis and exothermic 
combustion, but it is a process that requires fine tuning of 
operating conditions and feeding compositions. One parame¬ 
ter that is commonly introduced is the equivalent ratio (ER), the 
oxidizing agent fed (air in the present work) in relation to the 
oxidizing agent required for the stoichiometric combustion of 
the condensed matter (dry biomass in the present work). Using 
this definition, pyrolysis is a zero ER process, while combus¬ 
tion has an ER equal to or greater than 1. Gasification has an ER 
between 0 and 1 and its value is fundamental to the efficiency 
optimization of the gasification process. 

In this paper we present the results of a feasibility study of 
a biomass gasification plant with an innovative spouted bed 
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Notation 

A 

surface area (section area) [m 2 ] 

Cd 

drag coefficient 

D c 

column diameter [m] 

D, 

inlet diameter [m] 

D s 

spout diameter [m] 

d p 

particle diameter [m] 

ER 

equivalent ratio 

Fd 

particle drag force F D = C D |p g d 2 u 2 [N] 

G 

fluid mass flux in spouted bed referring to col¬ 
umn diameter [kg/m 2 s] 

9 

gravitational constant = 9.8166 [m/s 2 ] 

C if 

dto 

static bed height in spouted bed [m] 

oif 

<fl a 

bed (annulus) height in spouted bed [m] 

Ojf 

(/bC 

cone height in spouted bed [m] 

OJf 

dbj 

height of fountain in spouted bed [m] 

c if 

Jlfb 

height of free board in spouted bed [m] 

(flm 

maximum spoutable bed height [m] 

OLP 

dt m f 

minimum fluidization bed height [m] 

HHV 

higher heating value 

ICE 

internal combustion engine 

LHV 

lower heating value 

P 

pressure [Pa] 

Pi 

partial pressure of i in gas phase [Pa] 

AP f 

fluidization gas pressure drop [Pa] 

APmax 

maximum bed pressure drop [Pa] 

AP S 

stable spouting pressure drop [Pa] 

T 

temperature [k] 

U 

actual spouting velocity referring to column 
diameter [m/s] 

u t 

particle terminal settling velocity [m/s] 

Umf 

minimum fluidization velocity referring to col¬ 
umn diameter [m/s] 

Ums 

minimum spouting velocity referring to col¬ 
umn diameter [m/s] 

Ui 

gas velocity at the inlet of spouted beds [m/s] 

%,0 ffl a 

particle velocity on the axis at the bed surface 
of spouted beds 

Vsp 

solid particle volume [m 3 ] 

Greek letters 

tan y 

internal friction factor 

£ 

void fraction (gas phase fraction) 

£ mf 

bed void fraction at minimum fluidization con¬ 
dition 

£ 0 

bed void fraction of static bed 

£ ms 

bed void fraction at minimum spouting condi¬ 
tion in spouted bed 


spout voidage at the bed surface 

11 

yield - efficiency 

Q 

particle sphericity coefficient 

© 

base diverging angle [°] 

m 

gas viscosity [Pa s] 

Plop 

bulk density of solid particle [kg/m 3 ] 

Psp 

solid particle density [kg/m 3 ] 

Pg 

gas phase density [kg/m 3 ] 

Adimensional numbers 

Re 

Reynolds 

Re mf 

minimum fluidization Reynolds number Re m ^ = 

U m f d p pg 
p. 


Re t particle settling Reynolds number Re t = Utdppsp 
A r Archimedes Ar = dp ^ Psp Pg ^ pg9 

Li z 


reactor. The process was simulated through a thermodynamic 
approach using LIBPF, a C++ process modelling library. After 
a first simulation trial, a preliminary sizing of the plant was 
performed, revealing that the spouting regime of the gasifier 
imposed important constraints. A proper plant configuration 
was then developed and the main sensitivities were checked 
to study the optimal process conditions and performance. 

2. The spouted bed 

Gasification reactors, named gasifiers, are usually classified 
according to their solid handling characteristics. In order of 
increasing cost, complexity, efficiency and typical sizes, there 
are moving bed (downdraft and updraft), fluidized bed, bub¬ 
bling fluidized bed and entrained flow gasifiers (Higman and 
Burgt, 2003). 

In the present work we introduce an innovative application 
of spouted bed reactors. The spouted bed is a specific type of 
fluidized bed able to produce a good fluid-solid interaction. 
While in a standard fluidized bed, a fluid is passed through 
a uniform distributor plate to float the particles which move 
up and down in groups; in a spouted bed, instead, gas (or 
occasionally a liquid) is injected vertically upwards through 
a single central orifice into a bed of granular solid, expanding 
it. If the flow rate of the fluid and the pressure drop are suffi¬ 
cient and the bed depth is less than the “maximum spoutable 
bed depth”, the central jet breaks through the upper surface, 
resulting in a characteristic flow pattern known as spouting. 

The high fluid velocity causes a stream of solid particles to 
rise rapidly in a hollow central core within the bed, with high 
void fraction (~85%). The central core is called a “spout”, the 
peripheral annular region an “annulus” and the mushroom¬ 
shaped zone above the annulus a “fountain” (Fig. 1). The 
particles, having passed through the bed, fall back onto the 
annular space between the spout and the container wall form¬ 
ing a fountain, and then slowly travel downwards and, to some 
extent, inwards in a loosely packed bed. As the fluid travels 
upwards, it flares out into the annulus, which has a lower void 
fraction (40-50%) (Markowski and Kaminski, 1983; Bi, 2004; 
Anabtawi et al., 2005; Salam and Bhattacharya, 2006a, b). 

Spouted beds are the focus of renewed attention for physi¬ 
cal and chemical operations that require handling of solids for 
which fluidized beds present problems because of their great 
size, irregular texture, or sticky nature, as is the case of com¬ 
bustion of vegetable biomass or bituminous coal, pyrolysis of 
vegetable biomass, pyrolysis of waste plastics, and catalytic 
polymerizations (Markowski and Kaminski, 1983; Bilbao et al., 
1987; Rovero and Watkinson, 1990; Zahed and Epstein, 1992; 
Olazar et al., 1994, 2000; Aguado et al., 2000; Rasul, 2001; San 
Jose et al., 2004; Cui and Grace, 2007, 2008; Amutio et al., 2012; 
Arabiourrutia et al., 2012). 

In addition, spouted beds are of special relevance to 
biomass particle processing (drying, toasting, pyrolysis, gasi¬ 
fication and combustion). Some studies (Olazar et al., 1994) 
have demonstrated that spouted beds are quite good for A, 
B and D particle types of the well known Geldart’s classifica¬ 
tion (Geldart, 1973). Moreover, compared to standard fluidized 
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Fig. 1 - Spouted bed working scheme. 


beds, the spouted beds, and especially the dilute conical 
spouted beds, are suitable for avoiding defluidization due to 
bed agglomeration (Aguado et al., 2005). 

Gasification requires the development of technical pro¬ 
cesses that allow biomass particles to be spouted effectively 
with high gas-solid contact efficiency, high heat and mass 
transfer rates, as little flow of spouting gas as possible, high 
solid handling capacity, low particle entrainment and the 
shortest possible processing time. Gasification of solid fuels 
in a spouted bed, which has certain potential advantages 
over other fluid bed configurations, appears to be an under¬ 
exploited technique so far (Salam and Bhattacharya, 2006b). 

While gasification in a fixed bed can be considered as 
occurring in four different reaction zones (drying, pyrolysis, 
reduction and combustion) in a plug or sequential manner, in 
fluidized or spouted bed gasifiers separate reaction zones do 
not exist because of better mixing. All the processes of drying, 
pyrolysis, reduction and combustion can be regarded as taking 
place simultaneously throughout the reactor volume although 
the intensity of any particular process may vary according to 
its location (Salam and Bhattacharya, 2006b). 

3. Modelling approach 

In order to study high-temperature biomass-to-energy pro¬ 
cesses and simulate them, it is important to develop a proper 
modelling strategy. In this work, the thermodynamic approach 
was used and thus, concentrated parameter models, such as 
perfectly mixed volumes, were developed. 

Properly-customized, perfectly-mixed units can model 
gasifiers, heat exchangers, blowers, pumps, mixers, con¬ 
densers and internal combustion engines or fuel cells. Thus, 
the process can be modelled by suitably combining these 
units. 

Although spouted bed are characterized by short resi¬ 
dence time for the gas phase, we assume that thermodynamic 


equilibrium was reached in the gasifier, in view of the scopes 
of the present work, which are the study and the evaluation of 
the whole process and the survey of the main issues regarding 
it. 

3.1. Simulation tool 

The chosen simulation tool was LIBPF™ 1.0 (LIBrary for Pro¬ 
cess Flowsheeting), Greppi (2006), a process flow-sheeting and 
modelling tool arranged as a C++ library. 

The approach proposed in this chapter has been developed 
in the C++ programming language and in accordance with the 
implementation structure of the LIBPF library. The required 
low-level concepts are components and reactions. 

3.2. Components 

All main chemical species involved in the process from 
input to output were taken into account. Each chemical was 
characterized by its physical (density, composition, etc.) and 
thermodynamic (specific heat, enthalpy, vapour pressure, etc.) 
properties. As exposed in a previous work (Bernocco et al., 
2011 ), fifteen chemical species or components were taken into 
account in five conceptual component groups: dry and ash¬ 
free (daf) biomass, primary components (C, CH 4 , CO, CO 2 , 
H 2 , H 2 0, N 2 , NH 3 , 0 2 ), secondary components (methanol, 
ethylene), tertiary components (phenol, naphthalene), and 
inert/ash component (silica). Although methanol and ethyl¬ 
ene are unstable under gasification conditions, the choice of 
these secondary components is due to the representation of 
light hydrocarbons with a condensing temperature quite lower 
than that of water. Similarly, phenol and naphthalene were 
chosen to represent heavy tars with a condensing temperature 
higher than water. 

The last four component groups listed above are well- 
defined chemical species, so they were modelled using their 
standard LIBPF implementation. 

The biomass component requires a different modelling 
approach because it is not a well-defined chemical specie. The 
biomass model was developed as a chemical component with 
constant properties over all the calculations (flow sheet, unit 
operations, streams) similar to a well defined standard chem¬ 
ical species. Generic biomass is composed of moisture, ash 
and organic fraction (fixed carbon and volatile matter). Since 
moisture and ash can be taken into account separately, the 
biomass model involves only the organic fraction daf. 

3.3. Reactions 

The chemical species considered are involved in complex 
chemical transformation within a high-temperature process 
such as gasification. In this work three groups of chemical 
reactions were defined, and summarized in Table 1. 

The first two types are stoichiometrically defined for well- 
defined chemical species, while the pyrolysis reaction requires 
proper modelling of the biomass reactant. 

All high-temperature processing of biomasses, such as 
pyrolysis, gasification and combustion, develop hundreds of 
chemical species, and the stoichiometric relationship between 
products and biomass composition is complex to represent 
(Milne et al., 1998). Moreover the solid does not react directly 
with the gas phase except for solid carbon (char or coke). 
The solid biomass always undergoes a first step of pyrolysis 
or devolatilization (Souza-Santos, 2004). A systematic, generic 
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Table 1 - Chemical reactions considered. 


Combustions 

CH 4 + 20 2 —>■ C0 2 + 2H 2 0 

(i) 

H 2 + 0.5O 2 -> H 2 0 

( 2 ) 

CO 2 +0.502^ CO 2 

(3) 

C + O 2 —> CO 2 

(4) 

NH 3 + 0.75O 2 ^ 0.5N 2 + 1.5H 2 O 

(5) 

CH 3 OH + 1.50 2 -» C0 2 + 2H 2 0 

( 6 ) 

C 2 H 4 + 30 2 — ► 2C0 2 + 2H 2 0 

(7) 

Q 5 H 5 OH + 70 2 ^ 6C0 2 + 3H 2 0 

( 8 ) 

C 10 H 8 + 12O 2 ^ 10CO 2 +4H 2 0 

Equilibria 

(9) 

CH 4 + H 2 O ^ CO + 3H 2 

( 10 ) 

0.5N 2 + 1.5H 2 <=* NH 3 

( 11 ) 

co 2 + h 2 ^ CO + H 2 0 

( 12 ) 

2CO «=> C0 2 + C 

(13) 

C + O 2 ^ CO 2 

Reforming 

(14) 

CH 3 OH + H 2 0 -> CO 2 + 3 H 2 

(15) 

Q 5 H 5 OH + 5H 2 0 -> 6 CO + 8 H 2 

(16) 

C 10 H 8 + 10H 2 O -> 10CO + 14H 2 

(17) 

C 2 H 4 + 2H 2 0 -> 2CO + 4H 2 

Pyrolysis 

C c HhN n 0 0 l»h 2 H 2 + vq 2 O 2 + vqC + un 2 N 2 

(18) 

+ucoCO + uco 2 C 0 2 + uh 2 oH 2 0 + UCH 4 CH 4 

+l»nh 3 NH 3 

“5 ^c 2 h 4 C 2 H 4 + l»ch 3 ohCH 3 OH 

+l>c 6 h 5 ohC6H50H + l>c 10 h 8 CioH 8 

(19) 


stoichiometric approach to represent a biomass starting mate¬ 
rial and the repartition of pyrolysis products was developed in 
a previous work (Bernocco et al., 2011). The approach trans¬ 
forms the conventional representation of the starting organic 
matter of a biomass into a set of well-defined chemical species 
while respecting the atomic balances. The complex mixture 
of the possible gases and condensable products has been sim¬ 
plified (lumped) into 13 components and grouped as primary 
products (small molecules usually present in the pyrogas or 
syngas, also under thermodynamic equilibrium conditions), 
secondary products (light tars with a condensing temperature 
lower than that of water) and tertiary products (heavy tars with 
a condensing temperature higher than water). 

3.4. Unit operations 

Thermodynamic models generally present a system as a zero¬ 
dimensional point, with perfect homogeneity and immobility 
with respect to time. Thus the achievement of thermodynamic 
equilibrium was assumed. 

For these reasons the reactors, such as the gasifier, were 
modelled as ideal thermodynamic equilibrium reactors. 

The power consumption of the blowers was estimated by 
assuming isentropic compression efficiency equal to 70% and 
mechanical and electrical efficiencies equal to 90%. 

The cyclone was modelled as a perfect gas-solid separa¬ 
tor and sized using Stairmand’s procedure for high-efficiency 
cyclones (Sinnott et al., 1999). 

The syngas clean-up and cooling-condensing sections 
were modelled as two cooling units with an estimated pres¬ 
sure drop of 5500 Pa each. An outlet temperature of 130 °C was 
fixed for the first cooler (the clean-up unit): this involves the 
production of an high temperature usable heat, so that the 
clean-up unit was considered as a valuable heat cogenerator in 
the plant efficiency calculations. The syngas cooler-condenser 
was modelled as a cooling unit with flash vapour-liquid 



equilibrium with a fixed outlet temperature of 30 °C. This unit 
did not produce usable heat, and so was not included in the 
plant efficiency calculations. 

The power consumption of the pumps was estimated by 
assuming an electrical efficiency of 85% and a mechanical 
efficiency of 80%. 

The energy conversion section was represented as an inter¬ 
nal combustion engine, but modelled as an ideal combustor 
with outlet temperature fixed at 80 °C. The duty of this unit 
was multiplied by an electricity conversion efficiency of 30% 
in order to calculate the electrical power production. The 
remaining duty was multiplied by 75% efficiency in order to 
calculate the cogeneration heat power of the ICE (Internal 
Combustion Engine). These efficiency values are commonly 
accepted for small internal combustion engines, but they can 
also be extended to other technologies such as gas turbines. 
Bigger internal combustion engines for stationary combined 
heat and power generation can reach 48.7% of electrical power 
efficiency and 41.3% of cogenerating heat efficiency (for a 
9.5 MW Jenbacher engine). 

3.5. Key spouting parameters 

The available literature on spouted beds is mostly based 
on cylindrical apparatuses with diameters less than 0.3 m 
(Epstein and Grace, 2011). Studies of systems of lm diame¬ 
ter or larger are rare (Cui and Grace, 2008). So correlations and 
models mostly refer to small-scale vessels. 

The properties of spouted bed devices are well described 
by a few parameters (Becker, 1961), which are classifiable into 
three categories based on dependence. Fig. 2 should clarify 
the geometrical characteristics and the table in Appendix A 
summarises the major correlations for conical-based spouted 
beds. 
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A first parameter category can be considered as involv¬ 
ing parameters which depend only on the properties of the 
gas and particles and the process specifications. It groups 
the particle terminal settling velocity, hold up and minimum 
superficial velocity for fluidization. 

The terminal settling velocity of free-falling particles under 
gravity in a stagnant fluid is defined as the velocity at which 
the drag exerted by the surrounding fluid exactly counterbal¬ 
ances the net gravitational (weight minus buoyancy) force. 

The process specifications and average internal conditions 
give information about hold ups of solid (solid residence time) 
or about gas hour-space velocity (GHSV). This information 
usually comes from general kinetic approaches. 

The minimum velocity for fluidization is commonly 
defined as the superficial fluid velocity at which the extrap¬ 
olated relationship between the frictional pressure drop and 
the fluid velocity of a quiescent bed intersects the ordinate 
of equality between the drag and weight forces. The defini¬ 
tion is made unique by specifying that the quiescent bed is 
consolidated by a gradual reduction of agitation. 

A second category can be considered as involving parame¬ 
ters that also depend on the column geometry and its groups: 
inlet and column diameter constraints, maximum spoutable 
bed height, static and actual bed height and average spout 
diameter. 

The inlet diameter, column diameter and particle size are 
strictly related to spoutability and stability conditions deriving 
from experienced observation and experimental studies. 

There is a maximum bed depth which can still allow a sta¬ 
ble spout. The maximum spoutable bed height is an important 
parameter in the design of a spouted bed, because it is directly 
related to the amount of solid that can be processed. Above the 
maximum spoutable height, the bed enters bubbly and slug 
flow regimes. Transition to minimum fluidization occurs due 
to the gas diffusion in the packed bed. 

The bed depth can be defined as static or actual. The 
former comes directly from the solid bulk density and the col¬ 
umn geometry. The second can be estimated by using void 
fractions and spout and fountain sizes. However only exper¬ 
imental proof can correctly evaluate the actual bed height, 
by measuring the settled bed depth, observed after strongly 
agitating a bed and then gradually shutting off the gas flow 
(Becker, 1961). 

The spout diameter and shape at the bed level are only 
important at the gas inlet zone of the bed. Consequently, to 
model the gas and solid flow in spouted beds and for the sake 
of simplicity, it is useful to assume that the spout diameter 
is independent of its longitudinal position on the bed (San 
Jose et al., 2001). Usually, an average value of the diameter, 
estimated from semi-empirical or theoretical correlations, is 
taken. 

Finally, a third category can be considered as involving 
parameters that also depend on the bed depth, which is to 
say that it is related to the residence times of the gas or solid. 
This category groups the minimum and actual spouting veloc¬ 
ity, fountain size, maximum bed pressure drop, and stable bed 
pressure drop. 

The minimum spouting velocity is the lowest gas velocity 
required to make the fountain. The actual spouting velocity 
can be calculated from the gas feed rate, referring to the inlet 
section or to the column section. A slight reduction in the air 
velocity under this condition causes the spout to collapse and 
the pressure drop in the bed to rise suddenly. The minimum 
spouting velocity has been widely studied in the literature 


(Mathur and Epstein, 1974; Wu et al., 1987; Ye et al., 1992; Salam 
and Bhattacharya, 2006a). 

The fountain is characterized by a generally paraboloidal 
external shape. So it can be described by only two param¬ 
eters: the height and the base surface. The fountain height 
is a fundamental variable to evaluate the total height of 
the apparatus, while the base surface is important for the 
fountain-annulus interaction and for the solid recirculation 
(Wu et al., 1987; Ye et al., 1992) 

The maximum bed-pressure drop is the gas-pressure drop 
required to start up the spout and it is fundamental for sizing 
the gas-blowing system. 

The stable bed pressure drop is the gas pressure drop 
required for spouting under a stationary regime. 

In order to improve the spouting stability and gas flow rate 
flexibility, there are some works about the introduction of a 
draft tube that stabilizes the spout-annulus interface (Altzibar 
et al., 2005, 2009). In the present work, the draft tube was not 
considered in order to simplify the spouted bed design. 

Considerable information exists on the hydrodynamics of 
spouted beds and spout-fluidized beds under ambient condi¬ 
tions (Mathur and Epstein, 1974; Epstein and Grace, 2011). Only 
a few studies have been conducted on hot spouting from 20 to 
400°C (Wu et al., 1987) and from 20 to 880°C (Ye et al., 1992), 
applicable to combustion and gasification in conical-based 
cylindrical spouted beds. Some recent studies were done on 
high temperature spouting (Olazar et al., 2009; Makibar et al., 
2011 ), but they are strictly related to conical spouted beds. 

So there is considerable uncertainty that necessitates care 
and overdesign in sizing spouted bed gasifiers. Equations orig¬ 
inally developed for room-temperature conditions are often 
applied to high temperatures, with the assumption that these 
equations will give reasonable predictions if used with gas 
physical-property values appropriate to the actual operating 
conditions (Wu et al., 1987). 

For the sake of the applicability of the process study, the 
spouted bed sizing was important because it introduces the 
constraint of volumetric flow fed to the gasifier and so it 
requires the proposed recycling configuration. 

4. The reference process 

In the distributed energy production scenario the smallest 
plant capacity compatible with technological constraints is 
preferable. So, in choosing the plant size for our case study, the 
commercial availability of the biomass feedstock, the present 
knowledge of the spouted bed and the scope of energy pro¬ 
duction were considered. In our work, we assumed a nominal 
size of about 100 kW total thermal power, but the results can 
be generalized to smaller or greater sizes. This value should 
correspond to a biomass feeding rate of about 20 kg/h. Increas¬ 
ing the nominal plant capacity is possible, but not compatible 
with the aim of distributed energy production. On the other 
hand, reducing the plant size is not safe for the handling of 
commercially available feedstocks. It is possible to obtain a 
very low power-plant size using sawdust. This case involves 
such a small apparatus that the heat dispersion of the hot 
devices may undermine the feasibility of the unit operations 
from a thermal point of view. 

In a biomass gasification process, the CBL (Carbon 
Boundary Limit) delimits, on the basis of thermodynamic 
equilibrium, the region where solid carbon is all converted 
in gaseous products. The CBL strongly depends on ER, since 
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the latter influences the system composition by imposing the 
amount of oxidizing agent. For this reason, the optimization of 
the ER parameter is very important, as demonstrated by Prins’ 
exergetic analysis (Prins, 2005). In this work, after a few model 
test runs, we found that the best plant performance should be 
for ER varying between 30% and 50%. 

Another important parameter to be taken into account is 
the gasification temperature. Prins (2005) has demonstrated 
that about 927 °C must be reached for good reaction kinet¬ 
ics. Thus we adjusted the parameters in order to always reach 
930 °C in the spouted bed gasifier. 

The biomass feedstock considered in this study consisted 
of wood chips, whose nominal corn-position was taken from 
the Phyllis database (ECN, 2010) and was averaged as shown 
in Table 2. 

Commercially available wood chips are sold in four sizes, 
classified as P16, P45, P63 and P100 (European Standard 
EN14961-1). The energy consumption and the price of wood 
chips usually increase with decreasing size. For the sake 
of simplicity and as a compromise between solid dynamic 
requirements and commercial availability, wood chips of 
20 mm nominal size were chosen. This size can be classified 
as the average size of the P45 class of the European Standard. 

In order to evaluate the plant device sizing, the wood chips 
were further characterized by a particle density of 900 kg/m 3 
and a bulk density of 400kg/m 3 . 

The whole process is composed of five stages: air and 
biomass feeding, gasification, syngas clean-up and cooling, 
recycling, and electrical energy conversion as represented in 
Fig. 3. 

The biomass feeding system was assumed to be composed 
of an atmospheric hopper connected to a pressurized hopper, 
which fed a screw conveyor, through a rotary star valve. Using 
the model proposed by Dai and Grace (2008), the power con¬ 
sumption of the biomass feeding system was evaluated and it 
was so small as to be negligible. 

The air for the gasification was taken from the environ¬ 
ment and blown in the gasifying medium mixer using a proper 
blower. 

The gasifier was designed as a square-section spouted bed 
with a pyramidal base. The main parameters of this device 
were calculated using the correlations listed in Appendix A 
and are summarized in Table 3. 

The column side dimension, inlet diameter, base angle, 
base inlet side and cone height were designed according to 
correlations and best practice available in literature. The solid 
biomass residence time was estimated using the method pro¬ 
posed by Desappa et al. (2004) and the biomass content was 
consequently calculated. The inert (sand) mass fraction and 
content were fixed to a value which involves a sort of min¬ 
imum for the spouting velocity. The annulus voidage was 
assumed equal to the loosely packing voidage, which was 
experimentally measured using few material sample (wood 
chips and sand mixtures). The static bed depth was calculated 
using solid mass content and its bulk density and the geom¬ 
etry of the spouted bed. The maximum spoutable height, the 
average spout diameter, the minimum spouting velocity and 
the stable spouting pressure drop were calculated using the 
correlations listed in Appendix A. 

The syngas leaving the gasifier passes through the cyclone 
which produces a perfect ideal gas-solid separation by remov¬ 
ing ash and carbon. The resulting syngas passes to the 
clean-up section, represented as a simple cooler with an outlet 
temperature of 130 °C. 
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Fig. 3 - Process block diagram. 


After that the syngas is cooled further and the condensable 
component can be removed in the syngas cooler-condenser. 
We chose 30 °C as the outlet temperature of the syngas cooler- 
condenser as this value is easily reachable with either fresh 
water or with ambient-air cooling. 

The recycling section was developed in order to reach stable 
spouting gas velocity in the spouted bed. Different configu¬ 
rations were tested in our simulations. A first configuration 
without recycles was tested and demonstrated the need for a 
supplementary gas source to be fed to the spouted bed for a 
stable spouting regime. Adding more air was found to oxidize 
the produced syngas too much lowering its LHV. In a second 
configuration we tried to add low-pressure steam, but the heat 
cost for the steam production was too high for the plant power. 
A third configuration provided with flue-gas recirculation from 
the ICE was tested. Adding flue gases to the gasifier notably 
lowered the syngas LHV. This configuration could be exploited 
for a power generator that functions with a lean fuel, such as 
turbine gas. However, all of the tested configurations produced 
an outlet stream of syngas condensate. It was then proposed 
to recycle as much of the condensate as possible to reduce 
clean-up problems. 

The plant configuration described in the block diagram of 
Fig. 3 was found to be the best of all the tested configurations 
for gasification in a spouted bed reactor for energy production. 
This configuration consists of syngas condensate recycling 
and hot syngas recycling. The condensate recycling was 
designed firstly to limit waste-water production that exceeded 


Table 3 - Spouted bed gasifier parameters. 1 

Column side 

0.266 m 

Inlet diameter 

0.030 m 

Base angle 

60° 

Base inlet side 

0.040 m 

Cone height 

0.204 m 

Biomass residence time 

900 s 

Biomass content 

5.00 kg 

Inert particle dimension 

0 .001m 

Inert mass fraction 

0.75 

Inert content 

15.00 kg 

Annulus voidage 

0.506 

Static bed depth 

0.370 m 

Maximum spoutable height 

2.546 m 

Average spout diameter 

0.044 m 

Minimum spouting velocity 

0.464 m/s 

Stable spouting pressure drop 

3800 Pa 


the battery limits. But complete condensate recycling is dif¬ 
ficult to achieve, so the recycled condensate fraction is an 
important optimization parameter for the plant. 

The hot syngas recycling was added to adjust the volumet¬ 
ric flow fed to the gasifier. So the mass flow of the syngas 
recycling was constrained by the volumetric flow feeding the 
spouted bed acting on the stream split fraction after the 
cyclone. 

The condensate recycle and the syngas recycle contain 
combustible matter (methanol, phenol and naphthalene) that 
burns in the presence of hot air in the mixer before entering 
the spouted bed and in the inlet zone of the spout. So this 
recycling configuration should reduce the tar yield, but this 
desirable side effect could be confirmed only through a proper 
kinetic model of the gasifier. 

The purified and cooled syngas coming from the syngas 
cooler-condenser is passed to the power section, where it is 
burnt in the ICE to produce heat and electricity. 

5. Results 

A first modelling and simulation step was performed in order 
to obtain a preliminary stream table and a screenshot of the 
unit operation conditions. Our preliminary results were used 
to design the apparatus. The main result of this sizing was 
the determination of the volumetric flow constraint needed 
for the spouted bed gasifier. A value of 120 m 3 /h was found to 
be conservative for a stable spouting regime. 

Several important parameters can be simulated with the 
proposed model. The most important are the condensed 
recycle, and the air equivalent ratio, which directly influ¬ 
ence the efficiency of the process, so they represent process 
optimization parameters. As mentioned before, the gasifier 
was modelled as an adiabatic reactor, but introducing the 
conservative gasification temperature of 930 °C, ER and the 
condensate recycling were bound to that objective tempera¬ 
ture. 

Another parameter also influences the process behaviour: 
the quality of the feedstock. The dry biomass matter (organic 
matter and ash) composition is quite constant over a wide 
range of wood feedstocks, but the moisture can vary greatly. 
For this reason we studied the feedstock quality in terms of 
moisture content, varying this parameter from the ideal con¬ 
dition of 0wt% to the worst case of 35wt% moisture content 
in the biomass fed to the plant. 
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Biomass moisture fraction 

-0 - 0,05 - 0,1 - 0,15 - 0,2 - 0,25 - 0,3 0,35 



Fig. 4 - Mapping of the recycled condensate fraction (recycled condensate over total condensate from the syngas 
cooler-condenser) and air ER for different biomass moisture fractions. 


Many output variables can be monitored by the model, such 
as stream compositions and mass flows, pressures and tem¬ 
peratures, reaction conversions. For the scope of the process 
study, the focus was on the power and heat conversion. So the 
net power (electricity, cogeneration heat and total) yields refer¬ 
ring to the biomass feed rate and biomass LHV were monitored 
together with the syngas quality in terms of syngas LHV. 

Since the ER and condensate recycling are bound to the 
gasification temperature, it is possible to map a correspon¬ 
dence for the biomass moisture content as shown in Fig. 4. 
This map shows that for any given biomass moisture con¬ 
tent, increasing the condensate recycling requires a greater 
amount of air to be able to reach the gasification tempera¬ 
ture of 930 °C. This trend has two limits. The first consists of 
a thermal limit: when a high proportion of the condensate 
is recycled it becomes more difficult to reach the gasification 
temperature. The second consists of a mass balance limit: 
when too much condensate is recycled, the plant starts to 


accumulate a great amount of condensate until the steady 
state condition cannot be longer found. 

By increasing the biomass moisture content the exploitable 
region narrows in terms of condensate recycling. For the ideal 
0 wt% moisture content, it is possible to recycle all the con¬ 
densate with a limited quantity of air, otherwise, for the 
worst-case 35wt% moisture content, only about 45% of the 
total condensate coming from the syngas cooler-condenser 
can be recycled, while burning out about half of the biomass 
feed (air ER = 50%). 

The production of the plant in terms of net electric yield 
was studied and Fig. 5 shows the plant performance varia¬ 
tion against the air ER for various biomass moisture contents. 
Using the map in Fig. 4 the graph can be redrawn against the 
recycled condensate fraction. For the ideal 0% biomass mois¬ 
ture fraction a maximum net electrical yield of nearly 25% 
could be reached. Only for this ideal curve can a derivative 
sign change be noted at the ER value of 32%. The slight initial 


Biomass moisture fraction 

-0 - 0,05 - 0,1 - 0,15 - 0,2 - 0,25 - 0,3 - 0,35 



Fig. 5 - Net electrical energy yield vs. air ER for various biomass moisture fractions. 
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Biomass moisture fraction 



air ER 

Fig. 6 - Cogeneration yield vs. air ER for different biomass moisture fractions. 


yield increment in this case is due to the favourable effect of 
water on the reforming reactions. After this maximum and for 
all other curves, increasing the air flow, or, similarly, increas¬ 
ing the condensate recycling, reduces the yield because of 
the higher combustion of the available LHV of the biomass. 
Moreover, increasing the biomass moisture content lowers the 
yield. 

The heat cogeneration of the plant could be an interesting 
resource. As shown in Fig. 6, the cogeneration yield ranges are 
greater than corresponding electrical ones and it is possible to 
find a maximum cogeneration yield for all biomass moisture 
fractions. The first yield increase is due to the increased raw 
syngas flowing through the clean-up section (the first cogen¬ 
eration source). The succeeding yield decrease is due to the 
syngas LHV reduction and thus to the ICE total power reduc¬ 
tion (the second cogeneration source). Air preheating seems 
to have little effect on the overall heat balance. 

When the biomass moisture content is increased, the 
cogeneration yield decreases and the curves tend to flatten 
due to the heat required for vaporization. 

The total yield of the plant can be obtained by summing the 
electrical and heat yields since both refer to the biomass feed 
energy content. To present the results from a different point 
of view, the total yield shown in Fig. 7 was traced against the 
recycled condensate fraction for different biomass moisture 
fractions. Increasing the biomass moisture fraction causes the 


Biomass moisture fraction 



0 % 20 % 40 % 60 % 80 % 100 % 

Condensate recycle fraction 

Fig. 7 - Total (heat and power) yield of the plant vs. the 
recycled condensate fraction for different biomass moisture 
fractions. 


total yield to decrease as considerable energy is consumed in 
the moisture vaporization. 

For biomass moisture fractions ranging from 0wt% to 
10 wt% the total yield is low due to the cogeneration trend 
contribution. Above 10wt% biomass moisture content, the 
cogeneration yield trend flattens to produce the maximum in 
the total yield trend. 

6. Conclusions 

The feasibility of a biomass gasification plant was studied 
using thermodynamic models, and considering the applica¬ 
bility of a spouted bed reactor as innovative technology. After 
a first process simulation, a preliminary estimation of the 
sizes of the plant devices was made, highlighting an impor¬ 
tant issue regarding the applicability of a spouted bed gasifier. 
While, on the one hand, the optimal thermodynamic condi¬ 
tions require a certain amount of air and/or steam, on the 
other hand, the spouting regime requires a constrained range 
of gas flows. This issue justifies the complex recycling system 
(syngas recycling and syngas condensate recycling) studied. 

The introduction of the condensate recycling and the gasi¬ 
fication temperature constraint of 930 °C limited the plant 
parameters. Thus a plant functioning map was drawn for dif¬ 
ferent biomass moisture fractions and considering the air ER 
and the recycled condensate fraction. 

The plant power and heat yields were monitored to check 
the performance limits of our gasification plant configuration. 
An ideal value near to 25% of the net electrical efficiency was 
found to be reachable only for ideal dry biomasses. A wide 
range of values were found for the net electrical yield from 
the ideal maximum of 25% to a minimum of about 17% for the 
process conditions studied. The heat cogeneration yield was 
found to vary in the range 45-48% for all the process conditions 
studied. 

Further investigations are needed for plant sizing and func¬ 
tioning in order to develop a complete economic feasibility 
study. 
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Appendix A. 


Conical-based spouted bed design 


First category: dependence only 
from the properties of the gas and 
particles and process specification 


Second category: dependence also 
from column geometry 


Particle terminal settling velocity Ganser 
(1993) correlation (Ganser, 1993) 

Biomass particle hold ups 


Minimum fluidization velocity; Lippens 
and Mulder (1993) correlation 
(Coltters and Rivas, 2004) 


C D = ^(1 + 0.112[Re t K 1 K 2 ] a657 ) + 

l+(3305)/(Re t 2 K 1 K 2 ) ( Re t K l K 2 < 10 5 ) 

Few minutes for little particle, such as sawdust. 15-20 min for 
bigger particle, such as wood chips; The model proposed by 
Desappa et al. (2004) was used in the present work to 
estimate the residence time of wood chips 

Re^f = f( £ mf’ ®) + 9( £ mf’ ®)R e mf 


Minimum fluidization velocity; Wen and 
Yu correlation 

(Matsumura and Minowa, 2004) 

Becker condition 
(Mathur and Epstein, 1974) 

Spoutability, Chandnani and Epstein 
condition (Mathur and Epstein, 1974) 

Stability, Littman condition (Mathur and 
Epstein, 1974) 

Maximum spoutable bed height 

Malek and Lu equation (Mathur and 
Epstein, 1974) 

Maximum spoutable bed height McNab 
and Bridgwater (1977) equation 
(Ye et al., 1992) 


Re mf = [33.7 2 + 0.0408Ar]° 5 - 33.7 U mf = ^ 

(S)min ^-35 ^ or particles, 0.1 for thin particles 

< 25 30 Also valid for high temperature condition 

(Wu et al., 1987) 

ge > 3 4- 12 and 9€ a < 2£ m 


S 1 = ai05 (l) 075 (®f) 0 ' 4 $ : ( ^ in[Mgl) 


- D c (Dc\ 2/3 700 

- d p { Dj ) Ar 


(1 + 35.9 • 10 -6 Ar) ' - 1 


1 2 


Ki= (! + §0-°' 5 ) 1 K 2 = 
10 1.8148(-log 0) 0 5743 


Re mf = [C 2 + C 2 Ar]°' 5 - C i; Ci = 27.2, 
C 2 = 0.0408 recommended for high 
temperature conditions (Ye et al., 
1992) 


0 = 0.205-^3 (shape factor); D c = 

v p 

0.10^-0.23 m d p = 0.8 4 - 3.7 mm 

Suggested for high temperature 
conditions (Ye et al. 1992; Wu et al., 
1987). It over predicts Hm at room 
condition. The over prediction 
reduces increasing temperature. At 
~800°C, the over prediction is 
-30%. 
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Static bed height 


Average spout diameter 

Bridgwater and Mathur (1972) equation 

(Mathur and Epstein, 1974) 

Average spout diameter 
Wu et al. (1987) equation 
(Mathur and Epstein, 1974) 


Minimum spouting velocity 

Mathur and Gishler (1955) equation 

. , , , , (Mathur and Gishler, 1955) 

Third category: dependence also 

from bed depth 

Minimum spouting velocity 
Wu et al. (1987) equation 


Fountain height; Grace and Mathur (1978) 
theorical model (He et al., 1994) 

Fountain height; day theorical model; 
(Olazar et al., 2004) 

Maximum bed pressure drop 

Manurung (1964) equation 
(Mathur and Epstein, 1974) 


Stable bed pressure drop 
Manurung (1964) equation 
(Mathur and Epstein, 1974) 


Height deriving from solid 
residence time, solid bulk 
density in static condition and 
column geometry 

n _ 0.384 G°- 50 D°- 75 
U S ~ “025 

% 


D s = 


5.61G°- 433 d0- 58 V 133 

(P bp p g g) 0283 


U n 


= (2g%) a5 (|)(^) a33 («) C 


Suggested for high temperature 
conditions (Wu et al., 1987). It 
under exstimates Ds of ~30% (Ye 
et al., 1992) 

D c = 0.14 4- 1.1m i? = 60° ga = 

1.0 £= 1/6 d p = 

2.5 -T- 3.5 mm T = 30°C 


Ums — 


10 . 6 ( 2 g %) a5 (|) 105 (^) 0 ' 266 (^)“ 0 ' 095 («) 


0.256 


( £s p|%) 


1.46 (u 


Sp,0, dbg Psp 
2g Psp—Pg 


Suggested for high temperature 
conditions. It usually under 
predicts Ums by ~10% (Ye et al., 
1992; Wu et al., 1987) 


£=«•*(&-!) 


0.865 


) / tt tt \ 2.13 / \ —0.892 / \ - 

/ Pg u mf u pt \ I psp-pg \ I dp \ 

y Psp—Pg 9 D i J \ Pg ) y D i ) 


(S)-“ 


-iS = [^ + 08 °]- 34 - 4 4 ; 

QPbp^O 

S = 1 + ( 0 . 81 (tan y) 15 )/(f2)((D c d p )/(D i )) a78 ((D c )/(^ 0 ))’ 

f= __ 

d (same vo I. spherical particle) 


Suggested also for high 
temperature conditions (Wu et al., 
1987) 
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